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Coulomb repulsion is taken into account to derive the ther- 
modynamics of charged bosons in a random external poten- 
tial. A simple analytical form of the partition function is 
proposed for the case of non-overlapping localised states (i.e. 
a small amount of disorder). The density of localised bosons 
and the specific heat show a peculiar non-uniform tempera- 
ture dependence below the Bose-Einstein condensation tem- 
perature. The superfluid - Bose-glass phase diagram is dis- 
cussed. A new phase is predicted which is a Bose-glass at 
T — but a superfluid at finite T. 
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I. INTRODUCTION 

The charged Bose-gas (CBG), being an important ref- 
erence system of quantum statistics, has become of par- 
ticular physical interest in the field of high terpperature 
superconductivitytl. A long time ago SchafrothB demon- 
strated that an ideal gas of charged bosons exhibits the 
Meissner-Ochsenfeld effect below the ideal Bose-gas con- 
densation temperature. Later, the one-particle ij^cita- 
tion spectrum at T = was calculated by Foldy 
worked at zero temperature using the BogoliubovH 
proach. The Bogoliubov method leads to the result that 
the elementary excitations of the system have, for small 
momenta, energies characteristic of plasma oscillations 
which pass over smoothly for large momenta to the ener- 
gies characteristic of single particle excitations. Further 
investigations have been carried out at or near Tc, the 
transition temperature for the gas; Ihese have been con- 
cerned with the critical exponentsO and the change 
the transition temperature from that of the ideal gasB 
The RPA dielectric response function and screemng in a 
CBG have been studied in the high-density limitQiJ. The 
theory of the CBG beyond the lowest order Bogoliubov 
approximatiop-|Was discussed by Lee and FeenbergB and 
by Brucckncrll3, both of whom obtained the nextj-arder 
correction to the ground state energy. Woo and Matll cal- 
culated numerically the correctian to the Bogoliubov ex- 
citation spectrum. Alexandrovll3 found the critical mag- 
netic field Hc2{T) at which the CBG is condensed; the 
predicted temperature dependence of Hc2 was observed 
both in low-Tc and high- Tc cuprates, in which the coher- 
ence volume estimated from the heat capacity measure- 
ments is comparable with or even less than the unit cell 
volume. 

The doped Mott insulators are intrinsically disordered. 




Hence, the localisation of carriers in a random poten- 
tial plays a crucial role in their low-temperature ther- 
modynamics and transport properties. As an example, 
if a fraction of bipolarons is localised by disorder and 
the Coulomb repulsion in localised states is sufficiently 
strong, then the number of delocalised bosons is propor- 
tional to T while the boson-boson inelastic scattering rate 
is proportional to T^; this can explain both the linear 
temperature dependence of the in-plane [Sesistivity and 
the Hall density observed in the cuprateil. The picture 
of interacting bosons with a short-range interaction filling 
up all localised single-particle states in a random poten- 
tial and Bose-condensina into the first extended state is 
known in the literatureli3~tJ. To calculate the density of 
localised bosons one has to take into account the repul- 
sion between them. One cannot ignore the fact that the 
localisation length ^ generally varies with energy and di- 
verges at the mobility edge. One would expect that the 
number of hard core bosons in a localised state near the 
mobility edge diverges in a similar way to the localisation 
length. Therefore only a repulsive interaction can stop 
all particles condensing into the lowest localised state. 
Thus, as stressed by Fisher et aO, there is no sensi- 
ble non-interacting starting point in the Bose gas about 
which to perturb (in contrast to the Fermi gas). 

The scaling analysis p£ neutral bosons in random po- 
tentials by Fisher et a/t3 describes the Bose-glass - su- 
perfluid phase transition at zero temperature with the 
increasing density and (or) the hopping strength. How- 
ever the analysis is limiifid to the critical region near the 
transition and predictsEj the universal resistance R* at 
the Mott insulator to superconductor transition (in two 
dimensions) while dataHl on amorphous films suggest a 
wide range of R* . Simple analytically solvable models of 
interacting bosons in a rapdom potential might therefore 
be helpful. Lee and Gunnlij proposed a picture of neutral 
bosons in which the "true" extended Bose condensate co- 
exists with bosons in localised states. In this paper we 
develop a similar model for charged bosons interacting 
via Coulomb forces. 

The number of charged bosons in a single potential 
well is determined by the competition between the long- 
range Coulomb potential energy ~ 4e^/eo^ and the bind- 
ing energy Ec — E, where eg is the background dielectric 
constant and Ec is the mobility edge. If the Coulomb in- 
teraction is strong and the localisation length ^ diverges 
as ^ oc {Ec — E)^" with ly < 1, then each potential well 
cannot contain more than one bosonEI. Within this ap- 
proximation the localised bosons obey Fermi-Dirac statis- 
tics because Coulomb repulsion has the same effect as 
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the Pauli exclusion principle. In the extreme case of the 
hydrogen atom the average electron-nucleus distance is 
inversely proportional to the binding energy, i.e. v — I. 
The hydrogen negative ion exists but with rather low 
electron affinity so the doubly charged negative ion H^^ 
exists only as a resonant state. However, in general, the 
exponent v can be larger then unity and the Coulomb re- 
pulsion is not infinite, so the statistics of localised charged 
bosons is neither Bose nor Fermi. 

In this paper we study the parastatistics of partly lo- 
calised charged bosons in the superfluid phase. In this 
phase the chemical potential /i is exactly at the mobility 
edge, ji = Ec- We assume that the density of bosons n 
is not very high or the dimensionality is reduced so the 
localised states are not completely screened out (as in 
the case of the Mott transition in semiconductors). For a 
iD system that means that the dimensionless strength of 
the Coulomb repulsion should be about unity or larger, 
Ts = Ame^ / ^^{A'Knl'iY/^ > 1 where m is the boson mass 
and we take h — 1. At the same time the density should 
not be extremely low if delocalised bosons are in a super- 
fluid phase. A rather large magnitude of >> 1 is nec- 
essary for the Wigner crystallisation of charged carriers. 
Therefore, the above condition takes place practically in 
the whole relevant region of the carrier density. Moreover 
in reduced dimensions, relevant for cuprates, the localised 
states cannot be screened out even if < 1, and for 
rs > 1 the localisation length ^ is essentially unaffected 
by delocalised bosons or the Bogoliubov collective mode. 
Hence our picture suggests that some particles remain in 
localised states while others are in the "true" extended 
Bose condensate; for neutrial bosons this picture has been 
justified by Lee and GunnllS. 



the plasma frequencyu. Therefore the low-temperature 
thermodynamics is controlled by the excitation of the 
shallow localised states while the Bogoliubov collective 
modes can be ignored. In 3D their contribution is ex- 
ponentially small while in 2D their energy scales as 
and the specific heat as T^. Even in the case of a short 
range repulsion the sound modes yield an energy propor- 
tional to T''+^ and hence a specific heat which behayss 
like C cc T"^ (where d stands for the dimensionality )ll3. 
Therefore the contribution to thermodynamics from the 
delocalised bosons appears to be negligible at low tem- 
peratures compared with that from bosons localised in 
shallow potential wells for c? > 2. So in the following 
we calculate the partition function and specific heat of 
localised bosons only. 

When two or more charged bosons are in a single 
localised state of energy E there may be significant 
Coulomb energy and we try to take this into account as 
follows. The localisation length ^ is assumed to depend 
on E via 



(2) 



where v > Q. The Coulomb energy of p charged bosons 
confined within a radius ^ can be expected to be of order 



(3) 



Thus the total energy w oi p bosons in a localised state 
of energy E is taken to be 



w = pE + p{p — 1)k {—EY 



(4) 



II. PARTITION FUNCTION OF LOCALISED 
CHARGED BOSONS 

The Hamiltonian of charged bosons on an oppositely 
charged background (to ensure charge neutrality) in an 
external random field with the potential U (r) is given by 
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The Fourier component of the Coulomb potential for 
bosonic charge 2e is X^(k) = IGvre^/eofc^ in a ?>D sys- 
tem, and F(k) = 87re^/eofc in a 2D system with a three 
dimensional interaction. For a low amount of disorder a 
single particle spectrum consists of localised discrete lev- 
els below the bottom of the conduction band E^- At some 
finite temperature, Tc, bosons are condensed at E — Ec 
so that ji = Ec- If the Coulomb repulsion is strong one 
can expect that each localised state below E^ is occupied 
by one or a few bosons. The excitation spectrum of the 
delocalised charged superfluid has a gap of the order of 



where k > 0. Hence, as mentioned above, we see that the 
behaviour of charged bosons in localised states can be 
thought of as intermediate between Bose-Einstein statis- 
tics and Fermi-Dirac statistics. When k = we have an 
equally spaced set of levels, i.e. Bose-Einstein behaviour, 
whereas when k = oo we have Fermi-Dirac behaviour 
since the only levels with finite energy are p = and 
p — I, and thus an exclusion principle is enforced. When 
< K < oo we have the intermediate 'parastatistics' that 
the level spacing Sw — » oo as p — > oo. This behaviour is 
closely analagous to that observed in the Gojalomb block- 
ade model which applies to quantum dotalj. 

We are primarily interested in the properties of the su- 
perconducting phase, namely the phase in which a true 
(extended) Bose condensate is present. The only state in 
which the true condensate can form is the lowest energy 
delocalised state Ec- Hence, when there is a Bose con- 
densate, the chemical potential /i must be equal to zero 
if we choose Ec — 0- 

We take the total energy of a set of charged bosons in 
localised states to be the sum of the energies of the bosons 
in the individual potential wells. The partition function 
Z for such a system is then the product of the partition 
functions z for each of the wells, and the system free 
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energy F = ksTlnZ is simply the sum of the individual 
free energies fc^Tlnz. Hence it makes sense to study the 
partition function for one localised state of energy E on 
its own. The free energy of all localised bosons is then 
given by 



F = kBT f dENL{E)lnz{E), 

J —OO 



(5) 



where Nl{E) is the one-particle density of localised 
states below the mobility edge. 

There is an important reference case when the tem- 
perature dependence of the specific heat can be readily 
established; this occurs when v = \ and /i = 0. In that 
case the repulsive Coulomb energy scales as E, Eq.(4). 
By introducing E/ksT as a new variable in the integral, 
Eq.(5), and assuming that the density of states Nl is en- 
ergy independent in the region of order T/k below the 
mobility edge one obtains 
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In the general case f ^ 1 and (or) Nl not energy inde- 
pendent the specific heat has a non-linear temperature 
dependence. In this paper we study the case of low dop- 
ing when A^L(e) = it-lS^E — e) where ul is the number 
of impurity wells with only one localised level in each of 
them. 



III. CHARGED BOSONS IN A SINGLE 
LOCALISED STATE 

We focus on the properties of a single localised state e. 
The probability for the state to contain p bosons is pro- 
portional to e^^'-'^^P^'^ where is the chemical potential 
and P = I/UbT. (We shall retain the possibility of a 
non-zero chemical potential until it starts to complicate 
the equations.) We can re-express w — pfi as 
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■PIJ, = K {-ey {p - pof - K (-e) 
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where 



Po = 



2K(-e)'^ 



(7) 
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Fig 1 shows a graph of (w— p/z)/[k(— e)"^] as a function of 
p. The partition function ^(e) for such a single localised 
state is 



p=0 



p=0 



l3K(-eY {p-pof 



(9) 
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The partition function is thus completely determined by 
the dimensionless parameters po and ksT /[k{—€Y]. The 
mean occupancy {p) is 



d\nz{e) 
dfjL 



and, when /x = 0, the specific heat capacity c is 



c = 0- 



,9^ Inz(e) 
a/32 • 



(11) 
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Truncating the series at 100 terms the calculated values 
of these quantities are shown in Figs 2 and 3. We now at- 
tempt to understand these results in more detail, looking 
separately at each temperature range. 

1. fcijT <«■,(-£)'' 

At low temperatures the partition function is dom- 
inated by the term with p closest to pq, i.e. the 
value of p giving the lowest value oi w — p/i, and 
so the mean occupancy (p) is an integer and goes 
up in steps as po increases, as seen in Fig 2a. The 
changeover in dominance from one term to another 
occurs when po is a half-integer, at which point the 
two lowest energy states are degenerate. 

So long as one term dominates the partition func- 
tion, the specific heat capacity c will be close to 
zero. However when po is close to a half-integer we 
have a two level system and a corresponding Schot- 
tky anomaly in the specific heat capacity. This is 
seen in Fig 3a: c is zero when po is equal to a 
half-integer and rises to a maximum on either side, 
when the level separation is ~ fc^T. Hence, at 
fixed ksT /[n{—f'^)], the low temperature specific 
heat capacity (i.e. A;bT <C Pok(— e)'") is periodic 
in Po. We also note from Fig. 3 that the peak in 
the specific heat capacity rises to a maximum when 
Po is an integer. This is because the lowest energy 
level is the only non-degenerate level, all others be- 
ing in degenerate pairs; hence, at low temperature, 
we effectively have a two level system in which the 
upper level is a degenerate pair, thus resulting in 
a larger peak in c than occurs in an ordinary two 
level system. 

2. keT > K{-eY 

We can approximate the sum by an integral 
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(a) k{-€Y < kBT<plK{-eY 

In this case we can approximate the lower 
limit of the integral as —00, i.e. the partition 
function can be approximated by an untrun- 
cated gaussian, and is therefore approximately 
symmetrical about po. Hence, in this temper- 
ature range we have 



(14) 
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as is clearly seen in Figs 2a and 2b. 
After integration, the partition function be- 
comes 



from which we obtain that 



(15) 



(16) 



This result is another example of the equipar- 
tition theorem of classical statistical mechan- 
ics and arises simply because w is proportional 
to the square of a co-ordinate, namely p — pq. 
In Fig. 3b the effect of this result is seen in 
the region around ksT w k(— e)". 

(b) keT > K{~ey and kgT > plK{-eY 

In this case we can make the approximation 
that 
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and so the partition function becomes 



1^ 



(17) 



(18) 



As a consequence, the mean occupancy {p) be- 
comes proportional to the square root of tem- 
perature: 



{P) 



I ksT 
TTKl—eY 
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and this gives rise to the increase in (p) at 
large temperatures seen in Fig. 2. Once again, 
the specific heat capacity c ~ ^ks, as seen in 
Fig. 3. 



IV. BOSE-GLASS - SUPERFLUID TRANSITION 

Here we assume that all the localised states have the 
same value of e and derive the conditions for the Bose 
condensate to occur. The most startling conclusion is 
that, in particular circumstances, it is possible to take a 
system in which there is no Bose condensate and create 
one by raising the temperature. 

If the de-localised bosons are treated as being free par- 
ticles of spin-0 and mass m, then at temperature T and 
when ^ = 0, the number of de-localised bosons per unit 
cell of volume il is 



where 



knT, 



T 

Tco 



3.3125^ ( 1 



(20) 



If there are localised states in each unit cell and the 
mean occupancy of each state is (p) then the number of 
localised bosons per unit cell is 

If there are n bosons, in total, per unit cell, then the 
number of bosons in the Bose condensate is 

and the temperature of the phase transition Tc is found 
by solving 



■ riLip) 



Tea 



= 0. 



(21) 



However, as we have seen, {p) is in general a function of 
T and so solving this equation is not trivial. Nevertheless 
we can make some simple observations about the condi- 
tions necessary for a Bose condensate to occur. Fig. 4 
summarises the conclusions. 

When po is just above an integer value p_ , (jp) increases 
monotonically with temperature (as in Fig. 2b); we then 
have two possibilities: 

1. n < tt-lP- 

No Bose condensate occurs at any temperature. 

2. n> riLp^ 

In this case equation Eq.(|l]) has one solution, be- 
ing the temperature above which the Bose conden- 
sate disappears. 

When po is just below an integer value p+, on the other 
hand, (p) no longer increases monotonically with temper- 
ature (again, see Fig. 2b); rather, as T increases above 
zero (p) falls from p+ to po before rising back up again. 
Now we can distinguish three cases: 

1. n < riLPo 

No Bose condensate occurs at any temperature. 

2. n> ulp^ 

At T = the Bose condensate is present, but dis- 
appears as the temperature is raised. 
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3. UlPo < n < nLP+ This is more complicated. At 
T = the Bose condensate is absent. If 



kBTco <'■ 



then the Bose condensate wih be absent at all tem- 
peratures. If, on the other hand, 



then a Bose condensate will appear at a tempera- 
ture fcsT ~ k(— e)'" and disappear as the tempera- 
ture is raised further. 



V. CONCLUSION 

We have formulated the statistics of charged bosons 
partly localised by impurities. A simple form of the par- 
tition function is proposed by the use of a reasonable scal- 
ing of the Coulomb energy with the localisation length. A 
non-uniform dependence of the specific heat of the partly 
localised charged superfluid is found at low temperatures 
which strongly depends on the exponent v of the local- 
isation length. The Bose-glass - superfluid transition is 
analysed (Fig. 4) as a function of the n/riL ratio and the 
Coulomb interaction as characterised by the parameter 
Po- A new phase is found, which is a Bose-glass at T = 
but a superfluid at finite temperatures. 

We believe that our findings are relevant for doped 
high- Tc cuprates having-piany properties reminiscent of 
the charged Bose-liquidlll. Because the level of doping 
of these Mott-Hubbard insulators is high, one can expect 
Nl (E) to vary smoothly with energy E (rather than hav- 
ing spikes corresponding to discrete levels). In that case 
our approach leads to a non-linear temperature depen- 
dence of the specific heat at low temperatures ii v 1. 
From a preliminary analysis, the experimental observ^u 
tion of C oc T^~^ in superconducting La2-xSrxCu04t3 
appears to be consistent with the existence of bosons 
partly localised by disorder. 
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FIG. 1. Graph of {w — pfj,)/[K{—e)"] against p. 



FIG. 2. The mean occupancy (p) as a function of po and 
log^olkBT/lKi-e)"]}. (a) The full 3-D plot, (b) (p) versus 
log 10 {^sT'/[«;(— e)"]} for selected values of po- 



FIG. 3. The specific heat capacity c as a function of po 
and log^o {^BT/lK^-eY]}. (a) The fuU 3-D plot, (b) c versus 
logjQ {A:Br/[/«(— e)"]} for selected values of po: 3.0, 3.2, 3.3, 
3.4, 3.48 and 3.5. At low temperatures, c is approximately 
invariant under po — > po ± 1 and under po — > n — po where n 
is any integer greater than po. 

FIG. 4. Dependence of the possibility for Bose conden- 
sate formation on n/riL and po. In the shaded regions, no 

Bosc condensate can exist at T = but as the temperature 
is raised a condensate may form depending on the value of 

nlkBTco/H-ey]. 
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